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LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


The  contamination  of  spacecraft  surfaces  from  ambient  and  outgassing 
molecules  is  considered.  A model  is  formulated  for  the  molecules  desorbed 
from  the  spacecraft  surfaces  and  scattered  by  the  ambient  molecules  back  to 
the  emitting  and  adjacent  surfaces.  Equations  are  derived  for  the  backscatter 
and  far-field  cases  as  a function  of  surface  geometry,  mean  free  path,  and 
desorbing  beam  shape.  Expressions  are  given  for  point  to  point  and  flat  sur- 
face to  point  scattering.  The  scattered  molecular  irradiance  is  combined  with 
other  sources  to  obtain  the  total  molecular  irradiance.  A calculation  of  the 
contamination  during  launch  of  the  Lincoln  Experimental  Satellites  8 and  9 is 
given . 
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I. 


INTRODUCTION 


Spacecraft  contamination  due  to  the  adherence  of  molecules  on  external 
surfaces  can  in  some  instances  be  important  to  the  mission  success.  It  is 
well  known  that  the  buildup  of  molecular  layers  can  change  the  radiometric 
coefficients  of  a surface.  The  transmittance  and  reflectance  of  optical  sur- 
faces may  be  adversely  changed  in  terms  of  both  spectral  and  directional  charac- 
teristics. The  emissivity,  absorptivity,  and  reflectivity  of  nonoptical  sur- 
faces are  also  affected  by  surface  films.  These  films,  if  present,  must  be 
included  in  the  design  of  thermal  control  systems  and  sensors. 

In  this  Technical  Note  analytic  models  have  been  developed  for  estimating 

-2 

the  molecular  irradiance  (molecules  cm  ) on  a spacecraft  surface  due  to  mole- 
cules desorbed  (emitted)  from  adjacent  surfaces  and  scattered  by  ambient  mole- 
cules to  the  receiving  surface.  The  total  irradiance  on  a receiving  area 

(A  ) due  to  molecules  outgassing  from  an  emitting  area  (A  ) and  from  irradi- 
r e 

ating  ambient  molecules  in  the  free  stream  has  been  schematically  shown  in 
Fig.  1-1.  In  general,  the  irradiance  can  come  by  four  paths: 

(1)  Direct  irradiance:  molecules  are  desorbed  from  A^  and  go  directly 

to  A . In  this  case  the  surfaces  must  "see"  each  other  and  their 


r 


-+ 


normals  must  satisfy  the  condition  e • e >0. 

J e r 

(2)  Scattered  irradiance:  molecules  are  desorbed  from  A and  are  scat- 
e 

tered  by  another  surface  or  other  molecules  to  A^.  Models  for 
scattering  by  ambient  molecules  have  been  developed  in  this  Tech- 
nical Note. 
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Fig.  1-1.  Molecular  irradiance  paths. 
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(3)  Self  irradiance:  molecules  are  desorbed  by  A and  then  scattered 

j r 

back  to  A^.  Models  for  this  path  have  been  also  developed. 

(4)  Ambient  irradiance:  A is  bombarded  by  free  stream  molecules. 

r 

While  the  scattering  of  molecular  beams  by  molecules  and  surfaces  has 

been  intensively  investigated  both  theoretically  and  experimentally,  little 

work  has  been  done  on  molecular  exchanges  between  adjacent  surfaces  by  scat- 

1 2 

tering  with  intervening  molecules.  Scialdone  * has  developed  an  analytic 
model  for  the  self-contamination  of  a spherical  spacecraft  due  to  the  out- 
gassing  of  molecules  which  are  scattered  back  to  the  emitting  surface  (path  3) . 

In  the  present  analysis  this  path  has  been  considered  as  the  special  case  of 
the  emitting  and  receiving  areas  coinciding.  Scialdone  assumed  a simple  scat- 
tering model  for  the  desorbed  ambient  molecular  collisions.  A more  realistic, 
but  still  approximate,  isotropic  scattering  model  has  been  introduced  in  the 
present  analysis  and  it  is  shown  to  give  an  upper  bound  on  non-isotropic  scat- 
tering. Also,  it  is  shown  that  for  certain  cases  the  present  model  and  Scialdone's 

3 

model  give  identical  results.  Lee,  _et  al_ . reported  a Monte  Carlo  computer 

program  for  calculating  molecular  fluxes  within  vacuum  chambers.  The  program 

considered  only  scattering  by  the  chamber  walls  and  is  not  applicable  to  the 

problem  of  scattering  by  intervening  molecules.  The  problem  of  test  chamber 

2 

contamination  has  also  been  considered  by  Scialdone  using  his  backscatter 
model . 

In  Section  II  a general  equation  is  derived  for  the  molecular  irradiance 
assuming  single  isotropic  collisions.  The  resulting  volume  integral  is  evalu- 
ated in  Section  III  for  the  backscatter  case  (path  3)  and  is  compared  with 
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Scialdone's  model.  In  Section  IV  the  farfield  solution  is  obtained  for  the 
irradiance.  The  resulting  elementary  area  solution  is  integrated  to  get  the 
irradiance  at  a point  due  to  outgassing  from  an  adjacent  planar  area.  A sum- 
mary of  the  assumptions  and  formulas  is  given  in  Section  V.  The  resulting 
analytic  models  are  used  in  Section  VI  for  estimating  the  contamination  on  the 
Lincoln  Experimental  Satellites  (LES)  8 and  9 due  to  desorbed  molecules  from 
the  Transtage  launch  vehicle  during  boost  to  orbit.  Finally,  computational 
details  have  been  collected  in  the  Appendices. 

II.  IRRADIANCE  FROM  SCATTERED  MOLECULES 

As  shown  in  Fig.  1-1  the  molecular  irradiance  on  a surface  may  come  from 
ambient  molecules  directly  impacting  on  the  surface  or  from  adjacent  outgassing 
molecules  traveling  directly  to  the  surface  or,  thirdly,  from  adjacent  outgas- 
sing molecules  which  have  been  scattered  one  or  more  times  before  impinging  on 
the  surface.  In  any  particular  situation  several  of  these  irradiance  paths 
may  be  present.  The  non-scattered  components  of  the  total  irradiance  are  known 
once  the  surface  our gassing  rates,  ambient  molecular  fluxes,  and  the  geom- 
try  are  given.  The  scattered  component  is  more  complex  because  of  the  multitude 
of  paths  each  molecule  can  take,  the  details  of  the  collision  processes,  and 
geometry  factors  such  as  shadowing.  As  a first  approximation  it  will  be  assumed 
that  the  molecules  scattered  to  a surface  undergo  only  one  collision  with  an 
intervening  ambient  molecule  after  being  desorbed  from  an  adjacent  surface. 

This  model  is  applicable  to  the  low  density  regime  where  the  Knudsen  number 
(A/£)  is  large.  The  collision  processes  are  approximated  by  isotropic  scat- 
tering in  a reference  frame  attached  to  the  emitting  surface  (laboratory  frame). 
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This  assumption  simplifies  the  general  equation  in  Section  II. A;  as  discussed 
in  Section  II. B,  isotropic  scattering  also  gives  an  upper  bound  on  non-isotropic 
scattering  in  the  upstream  direction. 

A.  Isotropic  Scattering 

Consider  an  area  (A^)  emitting  molecules  into  the  solid  angle  du)  (Fig.  II-l)  . 

k 

By  definition  of  the  specific  intensity  (I)  , the  rate  of  molecules  leaving 

A in  direction  e is 
e h) 

= I A cos0  du).  (molecules  sec  S (1) 

e e 


The  total  emitted  molecular  flux  (N)  leaving  A^  is 


N = / 

a)  < 2tt 


-1  -2 

I cos0  do3.  (molecules  sec  cm  ) 
e e 


(2) 


For  a diffuse  or  Lambertian  surface  emitting  into  a conical  beam  defined  by 
the  semi-vertex  angle  0 , Eq.  (2)  becomes 


N = TTl  sin  0,  0 < 0 < tt/2. 


(3) 


Let  the  scattering  of  the  emitted  molecules  be  characterized  by  the  mean  free 
path  (A) . The  probability  of  a molecule  traveling  a distance  r^  without  being 

The  nomenclature  has  been  based  on  radiometric  symbols  with  molecules  used  in 
place  of  photons. ^ 

"k  k 

The  reason  for  considering  this  more  general  case  is  that  certain  integrals 
are  undefined  at  0=tt/2  and,  hence,  0 must  be  restricted. 
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Fig.  II-l.  Nomenclature  for  single  scattering  analysis 
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r 

scattered  is  by  definition  of  A given  by  exp(-  — The  number  of  molecules 

per  second  which  leave  A in  doo  and  arrive  at  distance  r is  therefore 

e e 


r 

0 

- I A cos0  do)  exp( t- ) . (4) 

e e A 

dr 

0 

The  fraction  of  these  molecules  scattered  in  distance  dr  is  — r — , hence,  the 

e A 

molecules  per  second  scattered  in  dr  is 

e 


r dr 
e e 

= I A cos0  do)  exp(-  -r-  ) — r — 
e e A A 


(5) 


The  scattered  molecules  in  general  leave  da)  and  travel  in  a new  direction. 
Assuming  isotropic  scattering  in  the  laboratory  frame f the  number  of  molecules 
per  second  entering  0)’  is 


r dr  f 

= I A cos©  doo  exp(-  ) ■ v-e  ~ . 
e e ^ X X 4tt 


(6) 


When  00 1 is  defined  to  contain  the  receiving  area  A^  (Fig.  II-l)  , the  molecular 
irradiance  (H)  on  A^ , assuming  no  further  scattering,  is 


A r dr  , _ ^ 

dH  = I cos0  da)  exp(-  ^ — . (molecules  sec  cm  ) (7) 

r 


Noting  that  du)  dr  = , where  dv  is  the  differential  scattering  volume, 

r 

0 

Eq.  (7)  becomes 


7 


dv 


(8) 


T A r 

jit  I e A , e x dv  f 

dH  = 7 — — cos0  exp( r-  ) 7 0) 

4tt  A e \ ' 2 

r Ar 


Equation  (8)  must  be  integrated  over  the  volume  (V)  defined  by  the  emitting 

71- 

beam  (0  < 0 < -^)  and  the  condition  that  0^  < tt/2  (disregarding  shadowing)  . The 
total  irradiance  becomes 


H 


HI  COs6e  eXp<-  T> 


(9) 


In  spite  of  the  simple  model  used  in  its  derivation,  Eq . (9)  is  difficult 
to  evaluate  for  general  orientations  between  the  emitting  and  receiving  areas. 
The  volume  V over  which  Eq.  (9)  is  integrated  has  been  indicated  in  Fig.  II-2. 
The  following  sections  evaluate  Eq . (9)  for  particular  cases  of  interest. 

B.  Nonisotropic  Scattering 

Appendix  A summarizes  briefly  the  theory  of  binary  molecular  collisions. 
For  an  elastic  collision  between  two  smooth  spheres  (adiabatic  collision)  it 
is  shown  that  the  scattering  of  both  spheres  is  isotropic  with  respect  to  a 
center  of  mass  coordinate  system.  Figure  II-3  shows  the  scattering  distri- 
bution of  a desorbed  molecule  after  collision  with  an  ambient  molecule  when 
both  molecules  have  been  approximated  by  perfectly  elastic,  smooth  spheres. 

It  can  be  seen  that  the  scattering  in  a reference  frame  attached  to  the  emit- 
ting area  (laboratory  frame)  is  nonisotropic  while  being  isotropic  in  a cen- 
ter of  mass  frame.  As  intuitively  expected.  Fig.  II-3  shows  that  more  de- 
sorbed molecules  are  scattered  downstream  in  the  direction  of  the  ambient 
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Fig.  II-2.  Integration  volumes. 
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MODEL  (centered  at  Vj)  OVER- 
BOUNDING UPSTREAM  SCATTERING 

Fig.  II-3.  Scattering  distribution  of  two  smooth  spheres. 
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motion  than  is  given  by  an  isotropic  distribution  in  the  laboratory  frame.  In 
Section  A the  fraction  of  molecules  scattered  into  a particular  solid  angle 
U)1  with  respect  to  a laboratory  reference  frame  has  been  approximated  by  the 
isotropic  distribution  The  distribution  will  give  an  upper  bound  on 

the  number  of  molecules  scattered  upstream  and  underestimates  the  number  scat- 
tered downstream;  hence,  the  irradiance  on  an  area  upstream  of  an  emitting 
area  is  upper  bounded  by  Eq.  (9).  The  nonisotropy  caused  by  non-spher ical 
molecules  is  difficult  to  estimate;  however,  this  effect  should  not  increase 
the  scattering  in  the  upstream  direction  so  that  Eq.  (9)  remains  an  upper  bound 
on  the  upstreams  irradiance. 

III.  BACKSCATTER  SOLUTIONS 

A.  Cylindrical  Beam  Backscattering 
Letting 

dv  * r ^ dr  doj, 
e e 


and 


A 

e 

A 


r 


= i. 


(10) 


Eq.  (9)  can  be  written  as 


H = 7 — / d(-r-)  exp(-  / cos9  da). 

4tt  j _ X p X \ \ e 
0 (Jd(r) 


(ID 
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Assuming  that  all  the  outgassing  molecules  are  confined  to  a narrow  cylindrical 


beam,  the  molecular  flux  N leaving  the  emitting  area  (A)  is 

N = I / cos0  doo,  (12) 

«o(r)  6 

where  Go(r)  is  the  solid  angle  enclosing  the  beam.  The  fraction  of  the  molecules 
isotropically  scattered  directly  back  to  the  emitting  surface  is  from  Eq.  (11) 

00 

S-i7/0d(f>  (13) 

The  scattering  solid  angle  go*  can  be  written  for  all  values  of  r as 


go  1 = 2tt(  1 - cosy)  , 


(14) 


where 


(15) 

(16) 

As  expected,  for  r large,  go?«-^  and  for  r small  go?~2tt. 

r 

Equation  (13)  is  nondimensionalized  by 


cosy  = 


<r2  + p2)1/2  ' 


TTp  = A. 
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X 


_r 

A’ 


e 


(17) 


Equation  (13)  can  be  written 


H 

N 2 


(18) 


where 


I,(e)  = / dx  exp(-x)  Tl  - — ~ „ /9 1 . (19) 

0 L (xz  + ez) 1/z  J 


* 

For  £ small,  it  can  be  shown  (Appendix  B)  that  I^(£)«£.  After  the  indi- 
cated substitutions,  Eq.  (18)  becomes 


H 

N 


* (' 


4ttA" 


1/2  A 

r)  , -o  « 1. 


(20) 


Scialdone  considered  molecules  radially  leaving  a spherical  spacecraft 
(radius  R)  and  being  scattered  back  toward  the  emitting  surface.  Molecules 
scattered  out  of  the  emitting  column  were  disregarded  under  the  assumption  that 


For  arbitrary  e the  exact  solution  is  a function  of  the  Struve  and  Neuman 
functions  of  the  first  kind.^  For  e small  I(£)  is  determined  by  the  integrand 
near  x=0;  hence,  bounding  exp(-x)  by  unity  and  integrating  gives  the  first  term 
in  the  expansion  of  the  exact  solution. 
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they  are  replaced  by  molecules  scattered  in  from  other  columns.  In  his  notation, 
the  returning  flux  NM  as  a function  of  the  number  of  original  molecules  per 


second  coming  from  the  surface  is 


(21) 


N. 


D 


The  departing  flux  N would  be 


4ttR 


— yielding  a flux  ratio  of 


= R 

N A 


(22) 


Using  Eq.  (20)  and  neglecting  scattering  from  other  beam  columns,  the  flux 
ratio  over  a spherical  spacecraft  also  becomes  y agreeing  with  Scialdone’s 
result.  A more  precise  analysis  taking  into  account  the  contribution  from 
other  beams  would  increase  the  flux  ratios  over  this  result.  This  extension 
is  difficult  as  will  become  apparent  from  considering  the  scattering  between 
two  surface  elements. 

B.  Conical  Beam  Backscatter ing 

Beginning  with  Eq.  (9)  with  = A^,  the  self-irradiance  of  a surface 


becomes 


H = Iff  cos6e  exP<~  f)  -^2 


V Ar 


(23) 
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The  integration  volume  V is  the  entire  emitting  cone.  Following  the  approach 
of  the  preceeding  section  let 


dv  = r sin  0 d0  d(f>  dr 
e e e 


(j0,=2tt 


1 - 


<rV>1/2 


TTy  = A COS  0 


Non-dimensionizing  as  before,  Eq.  (23)  becomes 


H = ttI  / d0  cos  0 sin0  l^(e,0  ) 
; « e e 2 e 


where 


I, (e,  6.)  = / dx  exp (- 


x)  |l 


. 2^  2 . .1/2 
(x  + e cos0  ) 
v e . 


By  arguments  similar  to  those  used  in  approximating  I^(e), 


1/2 

I0  (e,  0 ) cos  7 0 
2 e e 


Equation  (25)  becomes 

U 2 5/2 ..  ,ttA.  1/2 

(1  - cos  0)  K— ) 


Using  Eq.  (3) 


5/2 . 


H_2  (1-cos  ' 6)  , A ,1/2  A _ , 
sin2  6 (a’  V • 


(24) 


(25) 


(26) 

, 0 ) becomes 
e 

(27) 


(28) 


(29) 
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Two  limiting  values  for  Eq.  (29)  are  the  cylindrical  solution  0+0 


H 

N 


1/2 


0=0, 


and  the  hemispherical  solution 


(30) 


H 

N 


1/2 


0 = tt/2 


(31) 


A conical  beam  thus  produces  20%  less  backscattering  than  a cylindrical  beam 
with  the  same  emitted  molecular  flux. 
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IV.  FARFIELD  SOLUTIONS 

The  nomenclature  for  scattering  between  different  elementary  surfaces 
has  been  shown  in  Fig.  IV-1.  The  solid  angle  (aO  into  which  the  desorbed 
molecules  are  scattered  after  collision  with  an  ambient  molecule  is,  in 
general,  given  by  expressions  similar  to  Eqs.  (14)  and  (24).  A simplification 
is  possible  if  is  restricted  to  the  farfield.  In  this  limit 


and 


0) 


A cos0 
, r r 


1. 


Defining  the  nondimens ional  quantities 


r 


(32) 


(33) 


from  Fig.  IV-1,  cos  0^  can  be  written  as 


n S -*■  ->■  .X  -*■ 

cos  0 = — e • e + — e • e , 

r y s r y w r 


where 


2 2 2 ->  -> 
y = x + S 4-  2xS  e • e 
J a)  s 


(34) 


(35) 
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Fig.  IV- 1.  Orientation  between  two  elementary  surfaces. 
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Substituting  Eqs.  (32)  - (35)  in  Eq.  (9)  gives 


IA 

H - — ///  da)  dx  cosG^  exp(-x)  (—  eg  • er  + — e • e^).  (36) 

y s r y 


4ttX2  V 


The  volume  (V)  is  the  emitting  cone  for  which  cos  0^  >L  0.  Shadowing  by  other 

surfaces  has  been  disregarded,  but  could  be  included  if  necessary. 

A.  Arbitrary  Orientations  Between  Elementary  Surfaces 

When  the  volume  of  integration  is  the  entire  emitting  cone,  then 

Eq.  (36)  can  be  written  as 
IA 

e 


H 


4ttX 


//  do)  cos0e  (I3+I^)  , 


(37) 


where 


x3  - S 7 dx  , 


(38) 


i.  - « •«  7 <* x 

4 w r J 3 


(39) 


As  in  the  backscatter  solutions,  it  has  been  shown  in  Appendix  B that  bounding 

A 

I* 

exp  (-x)  by  unity  gives  a tight  bound  for  — r*  <<  1. 

Eqs.  (38)  and  (39)  can  be  integrated  to  give 


e *e 
s r 


l3~  ->  * 

S (1  + e • e ) 
v CO  s 


(40) 


e *e 

(0  r 


S(1  + eu-es) 


(41) 
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The  farfield  solution  becomes 


IA  (e*e+e*e) 

H = -£-  //  do  COS6e  — g-^—  r . (42) 

4ttAs  (1  + e *e  ) 

a)  s' 


As  shown  in  Fig.  IV-1  the  orientation  between  two  surfaces  can  be 

uniquely  described  by  the  distance  between  them  (s)  and  the  three  angles 

Yg 9 6g  and  0.  In  order  that  the  integration  be  over  the  entire  emitting  cone, 

the  emitting  surface  must  be  above  the  receiving  surface,  as  indicated  in 

Fig.  II-2,  and  the  emitting  cone  must  be  restricted  so  that  the  farfield 
A 

x* 

assumption  — — « 1 holds  over  the  entire  integration  range.  Under  these 
r 

limitations,  it  has  been  shown  in  Appendix  B that  starting  with  Eq.  (42)  the 
integration  over  the  entire  emitting  cone  gives 


H 

N 


- A 

A £ f 

4tt  As 


where 


cos  0 + cos  Y 
's 


1 + cos  0 

s 


(43) 


(44) 


This  result  is  independent  of  the  emitting  cone  angle  (0) . 

There  are  maximum  and  minimum  value  of  0 which  are  functions  of  0 and 

s 

Y . A consideration  of  the  possible  orientations  of  A with  respect  to  A 
fs  e r 

gives 


0 = y ”0?  0<y> 

min  ?s  s s — 's 


0 . =0  “ Y > 0>Y> 

min  s Ts  s — 's 
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0 =y+0,  0+v<7T, 

max  1 s s’  s ' s — 


0 


2tt  - (y  + 0 ),  0 + y 

1 « <5  ’ s ' S 


max  ' 1 s s 

A 

The  farfield  assumption  (—  « 1)  can  be 


>_  7T. 

written  using 


and 


so 


A 


r 


(r  ) 
r 


2 

min 


(r  ) = s sin (0  + 0 ) , 

r m s 


A 

r 

s2sin2(0  + 0 ) 
s 


« 1. 


(45) 


Figures  IV- 2 and  IV-3  show  the  function  f for  several  values  of  Yg* 

As  ys  approaches  0°,  f approaches  unity  and  H/N  becomes  independent  of  the 
surface  orientations  corresponding  to  the  emitting  surface  being  directly 
above  the  receiving  surface. 

B.  Plane  to  Coplanar  Point 

The  farfield  elementary  area  solution,  Eq.  (43),  can  be  applied 
to  find  the  molecular  irradiance  at  a point  due  to  desorption  from  an  adjacent 
coplanar  area.  Figure  IV-4  shows  the  nomenclature  for  a flat,  rectangular 
surface  (A)  emitting  molecules  which  are  scattered  to  a single  coplanar  point. 
In  this  case 

v = 0 = tt/2 

» s s 

“0=0, 

and 
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10-6  -16228 


Fig.  IV-3.  f as  a function  of  0 and  0 (y  = 45°,  30°). 

s s 
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Fig.  IV-4.  Scattering  from  a plane  to  a coplanar  point. 
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f = 1 


(46) 


Eq.  (43)  gives 


H = rr  dA 

N 4ttA  * { s' 

A 

Eq,  (47)  can  be  written  as 


H - fj  (y)  K(o,S), 


where 


a 1 

K(a,B)  = // 

1-1 


dx  dy 


, 2 . a2  2.1/2 
(x  + B y ) 


a 


r - b 
6 ■ i 


The  integration  of  Eq.  (49)  gives 

««.B)  = at. 

-B  + (a  + B ) 7 - 

. B + (1+  B2)1/2 
L-b  + d+  B2)1/2_ 


+23&n 


[a  + (a2+  e2)1/2~ 
|_1  + (1+B2)1^2 


Note  that  in  the  limit  of  OtKL  and  3“K),  Eq.  (50)  gives 
2B(a-l)  or  H = fjffx  where  A=2b  Figure  IV-4 

K(a,B). 


(47) 

(48) 


(49) 


(50) 

as  expected,  K(a,3) 
shows  the  function 
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V.  SUMMARY 

Basic  assumptions:  (1)  Lambertian  emitting  surface 

(2)  Single  collision  with  ambient  molecules 

(3)  Isotropic  scattering  in  laboratory  frame 


Backscattering : ~ = T c-°- — ( ~~9)  , 0 < 0 < ir/2 

N 5 sin20  ttA2  ~ ” 


I = A (JL_) 

N 5 4ttA2 


1/2 


, 0 =tt/2  (hemispheric  beam) 


-=  (— ) 
N \ ,2; 
4ttA 


1/2 


, &=0  (cylindrical  beam) 


Constraints : 


Fo)  « l 

A 


Backscattering  to  a sphere 
of  radius  R: 


H = R 
N A 


cylindrical  beam 


A = 4 Y hemispherical  beam 
N DA 


Scattering  from  adjacent  surfaces 
Farf ield  solution : 


H f e 

N = 4tt  ~Xs  ^n^ePen<^ent  °f  emitting 
cone  angle 


f = 


cos0  + cosy 


1 + COS0 


Constraints:  (1)  Receiving  area  "sees"  entire  emitting  cone 

A 


(2) 


(rr}  min 


« 1 
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O)  e . > e > 

min  — — 


max 


Scattering  from  a plane  to  a coplanar  point 


H 

N 


1_ 

4tt 


j K(o,e> 
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VI. 


EXAMPLE  CALCULATION  - LES-8/9 


During  the  launch  of  LES-8/9,  the  satellite  surfaces  will  be  exposed 
to  molecular  fluxes  from  several  sources.  Potential  contamination  of  the 
satellites  will  occur  (1)  during  enclosure  in  the  payload  fairing  where  the 
main  molecular  irradiance  source  will  be  the  inside  fairing  surface,  (2) 
between  payload  fairing  jettison  and  satellite  deployment  during  which  the 
main  molecular  irradiance  source  will  be  the  launch  vehicle  (Transtage),  and 

(3)  after  satellite  deployment.  This  analysis  estimates  the  molecular  irradiance 
between  payload  fairing  jettison  and  satellite  deployment  (a  period  of  about 
24,240  seconds)  with  the  goal  of  assuring  that  the  resulting  contamination 
will  not  degrade  the  LES-8/9  operational  performance. 

The  assumptions  used  in  the  analysis  and  discussed  in  detail  later  are: 

(1)  To  simplify  geometry  the  satellite  surfaces  are  characterized  by 
a "contamination  envelope".  Contamination  of  the  envelope  is 
greater  than  the  contamination  of  the  satellite  surfaces. 

(2)  All  irradiant  molecules  condense  on  the  contamination  envelope, 
i.e.,  none  bounce  off  or  are  re-emitted  (cold  surface  assumption). 

(3)  Surfaces  not  facing  each  other  may  exchange  molecules  by  scattering 
from  other  surfaces  or  ambient  molecules. 

(4)  The  scattering  from  ambient  molecules  is  given  by  the  model  in 
the  preceeding  sections. 

(5)  The  epoxy  paint  on  the  surface  of  the  Transtage  tanks  is  the  major 
source  of  contaminating  molecules. 


28 


(6)  The  rate  of  mass  fraction  loss  of  volatile  condensible  materials 
(VCM)  is  0.1%  per  24  hours  and  assumed  to  be  constant  with  time. 


(7)  Contamination  is  measured  by  the  cumulated  surface  density 

-2 

deposited  on  the  contamination  envelope  (gm  cm  ) . 

The  profile  of  the  launch  vehicle  (Transtage)  has  been  shown  in 

6 

Fig.  VI-1.  The  painted  surface  covering  the  tanks  has  been  indicated.  For  a 

_3 

paint  thickness  of  five  mils  with  a density  of  1.2  gm  cm  , the  amount  of  paint 
which  can  potentially  contaminate  is  about  six  pounds;  however,  only  a small 
fraction  of  this  will  actually  be  scattered  to  the  satellite  surfaces. 

Figure  VI-2  shows  details  of  the  payload  stack.  For  this  launch  the  payload 
consists  of  the  LES-8/9  satellites  and  two  SOLRAD  satellites.  Only  contamination 
of  LES-8/9  has  been  considered.  In  order  to  simplify  the  geometry  and  over- 
bound the  irradiance,  LES-8/9  has  been  replaced  by  a Contamination  envelope" 
with  the  shape  of  a cylinder  as  shown  in  Fig.  VI-2.  Since  it  is  closer  to  the 

Transtage  the  irradiance  on  this  envelope  by  molecules  originating  from  the 
sides  and  bottom  of  the  Transtage  over-bounds  the  irradiance  on  any  of  the 
LES-8/9  surfaces  from  the  same  sources. 

7 8 

Figure  VI-3  shows  some  details  of  the  transfer  trajectory.  ’ The  Transtage 
and  payload  are  injected  (482  seconds)  into  an  elliptic  park  orbit  of  about 

80  x 95  nmi  . at  the  end  of  the  State  II  depletion  burn.  The  Transtage  is 
ignited  near  parking  orbit  apogee  (4439  seconds)  and  produces  a Hohmann 
transfer  orbit  having  a synchronous  apogee.  The  Transtage  is  ignited  for  the 
second  time  (23,705  seconds)  to  obtain  a final  circular  synchronous  orbit  with 
an  inclination  of  22.84  degrees.  The  LES-8/9  satellites  are  deployed  (24,555 
seconds)  followed  by  the  SOLRAD  satellites  (25,330  seconds). 
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18-6-16230 


Fig.  VI-1.  Transtage  profile. 


30 


CONTAMINATION  ENVELOPE  |l8-6-1623~T1 


Fig.  VI-2.  Payload  stack  details. 
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BURN 


Fig.  VI-3. 


LES-8/9  transfer  trajectory. 
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Based  on  discussions  of  the  possible  contaminating  effluxes  from  the 
upper  stack,  it  has  been  assumed  that  the  painted  surface  of  the  Transtage  is 

•k 

likely  to  be  the  largest  source  of  VCM.  Since  this  surface  is  close  to  the 
LES-8/9  surfaces  and  is  heated  during  launch,  it  is  likely  to  be  a "dirty" 
surface  which  will  contaminate  to  a certain  extent  the  "clean"  spacecraft 
surfaces.  From  Fig.  VI-2  it  is  seen  that  the  dominant  path  of  molecules 
originating  from  surface  does  not  "see"  the  painted  Transtage  side.  Since 
the  Transtage  is  downstream  of  the  contamination  envelope,  the  theory  in  the 
preceeding  sections  will  over-bound  the  LES-8/9  contamination. 

For  the  altitudes  of  interest,  scattering  of  the  Transtage  molecules 
to  the  contamination  envelope  is  in  a free-molecular  regime,  i.e.,  there  is 
little  interaction  between  the  desorbed  Transtage  molecules  and  ambient 
molecules.  The  ratio  of  mean  free  path  to  a characteristic  dimension  (Knudsen 
number)  is  large.  For  example,  at  payload  fairing  jettison  X  *  **  13  meters  and 
for  l = 1 meter,  = 13.  At  higher  altitudes  the  Knudsen  number  is  further 
increased. 

The  relative  motion  between  the  launch  vehicle  and  ambient  molecules 

is  accounted  for  by  the  relative  mean  free  path.  The  relative  mean  free  path 

is  a function  of  the  angle  between  the  launch  vehicle  velocity  and  the  emitted 

molecule  direction. ^ When  this  angle  is  near  90  degrees,  which  is  true  for 

most  of  the  molecules  emitted  by  the  Transtage  tank  surface  in  parking  orbit, 

X ^ « X , . Hence,  all  calculations  have  assumed  ambient  mean  free  paths, 
rel  amb 

If  the  molecules  were  emitted  into  the  ambient  free  stream  then  A ^ < A 

* 

Jack  Lynch  and  Ed  Murphy,  private  communications,  April  1974. 
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Assuming  a Transtage  surface  temperature  of  100°C,  the  rate  of  VCM 
9 

mass  fraction  loss  (Y)  is. 


Y = 0.1%/24  hours  = 1.1574  x 10  ^ sec 


(51) 


The  molecular  efflux  (N)  of  VCM  from  the  Transtage  surface  is 


N = 


N 6pY 
o 

W 


(molecules  cm  ^ sec  , 


(52) 


where 

Nq  = Avogadro’s  number, 

6 = paint  thickness 

p = paint  density 
W = VCM  molecular  weight. 

The  contamination  is  measured  by  the  rate  of  mass  deposited  on  the  contamination 
envelope.  Assuming  all  the  molecules  condense  or  are  absorbed  on  the  envelope 
surface,  the  contamination  rate  (C)  is 
x _ HW  ( -2  -1. 

C = — (gm  cm  sec  ).  (53) 

o 

For  simplicity  the  curved  Transtage  tank  has  been  replaced  with  a flat 
rectangle  which  is  coplanar  to  the  irradiance  point.  The  area  of  this 
rectangle  has  been  taken  to  be  one  quarter  of  the  total  Transtage  tank  area. 

This  assumption  over-bounds  the  number  of  molecules  scattered  to  a particular 
envelope  point  since  the  Transtage  curvature  would  reduce  the  irradiance  below 
the  calculated  flux.  Since  the  effect  of  distance  on  the  irradiance  goes  as 
1/s,  Transtage  surfaces  closest  to  the  contamination  envelope  dominate  the 
irradiance. 
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Using  Eqs.  (48)  and  (52)  in  Eq.  (53)  gives, 
c - |{p  (f)  K(a,6)  (54) 

Note  that  C is  independent  of  the  molecular  weight  of  the  VCM.  The  cumulated 
contamination  is  the  time  integral  of  Eq.  (54), 
t . 

C(t)  = / C dt.  (55) 

o 

In  Eq.  (54)  all  the  parameters  have  been  assumed  independent  of  time 
except  X which  varies  according  to  the  trajectory.  Figure  IV-4  shows  the 
variation  of  X with  time  for  the  LES-8/9  launch.  The  maximum  envelope 
contamination  takes  place  closest  to  the  Transtage  (Fig.  VI-2)  and  when  X is 
minimum,  i.e.,  at  payload  fairing  jettison.  Table  1 lists  the  parameters  used 
to  compute  C and  C. 


TABLE  1 

PARAMETERS  FOR  LES-8/9  CONTAMINATION  COMPUTATIONS 

Transtage  Surface  Paint: 

6 = 5 mils 

-3 

p = 1.2  gm  cm 

Y = 0.1%  per  24  hours 

Geometry  of  Transtage/ Envelope : 


£ = 15  inches 

a = 5.96 

L = 89.7 

8 = 2.0 

b = 30.0 

K(a,8)=  6.365 
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ALTITUDE  (km) 


TIME  FROM  IGNITION  (sec) 


Fig.  VI-4.  LES-8/9  altitude  and  mean  free  path  histories. 
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MEAN  FREE  PATH  (m) 


With  the  above  values,  the  maximum  contamination  rate  immediately 
after  payload  fairing  jettision  is 

A 5x10  ^ x 2.54  x 1.2  x 1.1574x10  ® 0 ^ 

C = x 2.94x10  x 6.365 

4tt 

C = 2.638x10  cm  ^ sec  ^ . (56) 

Figure  IV-5  shows  a plot  of  the  contamination  rate  and  the  cumulated 
contamination  as  a function  of  time.  Because  of  the  sharp  increase  in  X 
with  time  (Fig.  IV-4),the  cumulated  contamination  remains  constant  after 
about  1400  seconds  from  ignition. 

The  conclusions  reached  by  this  analysis  are: 

(1)  The  maximum  contamination  rate  from  Transtage  surface  occurs 
immediately  after  payload  fairing  jettison  (305  seconds,  133  km  altitude). 

(2)  An  upper  bound  on  the  maximum  cumulated  contamination  is 
5.53  x 10-10  gms  cm-2  and  is  reached  at  about  1400  seconds  after  ignition. 

(3)  The  contamination  collected  between  payload  fairing  jettison  and 
satellite  deployment  is  not  expected  to  affect  the  operational  satellite 
performance . * 

^ ^ g _2  g 

A molecular  monolayer  is  about  10  gm  cm  (1^0  - 1.587x10  gm  cm  ). 
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CONTAMINATION  RATE  (gm  cm^sec-1) 


Fig.  VI-5.  LES-8/9  contamination  history. 
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CUMULATIVE  CONTAMINATION  (gm  cm‘2) 


APPENDIX  A 


BINARY  MOLECULAR  COLLISIONS 

Basic  Equations:  Consider  a molecule  (mass  m^)  desorbed  from  a 

surface  at  time  t=0  and  which  travels  in  a straight  line  with  respect  to  the 
surface*.  At  t=t  there  is  a collision  with  an  ambient  molecule  (mass  m,) 

c i 

resulting,  in  general,  in  a new  velocity  (v^)  and  a change  in  direction  (x^)* 
Figure  A-l  shows  this  scattering  process  in  a frame  of  reference  attached  to 
the  emitting  surface;  this  reference  frame  is  referred  to  as  the  laboratory 
frame  throughout  this  report. 

The  dynamics  of  colliding  molecules  are  described  in  part  by  Newton's 

10 

laws  of  motion  applied  to  collections  of  particles.  Summarizing  these  laws 
the  velocity  of  the  center  of  mass  (vcm)  is  as  if  the  total  system  mass  (M) 
were  concentrated  at  the  center  of  mass  (CM)  and  acted  upon  by  the  sum  of  the 
forces  (F)  external  to  the  system. 

F = Mv 

cm 


where 


F = E 

k=i,j 


Fk’ 


M = E M.  , 
k=i,j 


v = E 
Cm  k=i,j 


Vjj 

M 


(A-l) 


For  a system  of  two  colliding  particles  (m^,m^)  the  collision  forces  are 
internal  to  the  system  (F  = 0),  hence 


The  molecule  has  the  surface  velocity  plus  a velocity  with  respect  to  the 
surface  (v . ) which  for  this  analysis  has  been  assumed  to  be  distributed 
isotropically,  i.e.,  the  surface  is  Lambertian. 
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Fig.  A-l.  Geometry  of  desorbed/ambient  molecular  collisions. 
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Vk 


v = E 
cm  , 4 . M 
k=i,j 


= E 


M 


= constant. 


(A-2) 


k=i,j 

The  velocity  of  the  molecules  relative  to  the  CM  is, 

V,  = v - v , 
k k cm 

V'  = v."  - "v  , k=i,j 
k k cm  ,J 


(A- 3) 


By  definition  of  the  CM,  the  total  system  momentum  relative  to  the  CM  is  zero 
even  when  F ^ 0. 


E hlV  = E vlV'  = 0 
k=i,j  k R k=i,j 


(A-4) 


The  total  kinetic  energy  (T)  of  a system  can  be  written  as  the  kinetic 
energy  due  to  the  total  mass  moving  at  v plus  the  kinetic  energy  due  to 
motion  relative  to  the  CM. 


T = \ Mv2  + E \ mkVk 
2 Cm  k=i,j  2 


(A-5) 


When  F = 0 and  v = constant,  T will  be  decreased  during  the  collision  by 
cm 

internal  forces  moving  through  displacements  relative  to  the  CM.  The  total 
kinetic  energy  after  collision  (T^)  differs  from  the  intial  kinetic  energy 
by  the  energy  which  has  been  absorbed  by  the  internal  degrees  of  freedom  of 
the  molecules. 


T'  = i Mv2  + E im.v;2  = T-  AE 
2 cm  , , ,2  K k 


(A- 6) 


k=if  j 

When  AE=0,  the  collision  is  termed  adiabatic  and  corresponds  to  no  changes  in 

the  electronic,  vibrational  and  rotational  states  of  the  molecules. 

The  total  angular  momentum  of  a system  of  particles  relative  to  the  CM 

(H  ) is  related  to  the  moment  of  the  external  forces  about  the  center  of  mass 
cm  

cm 
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where 


M = H 
CM  CM  5 


m~  - z *k x Fk  • 


CM 


K-i.j 


H~-  z \ x Vk 


CM 


k=i,  j 


\ = r’  " r 


k CM  ’ 


- v Vk 


(A- 7) 


For  Fk=0, 


z \ X Vk  ■ . ( ®k  * Vk  ■ 0 

k=i,j  k=i, j 


(A- 8) 


Equations  (A-2),  (A-4),  (A-6)  and  (A-8)  are  not  sufficient  to  uniquely 
determine  the  resulting  velocity  (V^)  and  scattering  angle  (x^)  of  the  collision. 
A model  for  the  interaction  during  collision  must  be  added. 

Scattering  in  the  Center-of-Mass  System:  Figure  (A-2)  shows  the 

collision  geometry  in  the  CM  system.  Equations  (A-4)  give 


Vit  + Vjt  ■ Vlt  + ”jvit  - 0 ’ 

(A-9) 

When  the  interaction  potential  has  sperical  symmetry,  the  molecules  will  remain 


mV  + m.V , = M,V;  + m.V;  = 0 

i in  j jn  i in  j jn 


in  the  plane  defined  by  their  initial  velocities,  i.e.,  and  v ^ . 
The  smooth  sphere  model  gives, 
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b 


Fig.  A-2.  Collisions  of  hard  spheres  in  the  center  of  mass  system. 
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vkt 

= \t  ’ 

k=i,j 

To  solve 

for  v; 

and  V; 

let 

in 

jn 

- v;  = 

e(V . ■ 

jn 

in 

v in 

where 

0 < 

e <_  1 

• 

(A-10) 


The  coefficient  of  restitution  (e)  measures  the  degree  that  the  collision 
departs  from  an  adiabatic  one.  For  e=0,  the  collision  is  inelastic  (non- 
adiabatic)  and  for  e=l,  the  collision  is  perfectly  elastic  (adiabatic). 
Solving  Eqs.  (A-9)  and  (A-ll)  gives, 


from, 


V"  = -eV.  , k=i,j  . 
kn  kn  9 ,J 


(A-12) 


. th 


The  scattering  angle  of  the  i molecule  in  the  CM  system  is  found 


where 


V • V 
i i 

cos  X = — — > 

i i 


V,  = V e + V.  e,  , 
i It  t in  in 

V'  = V'  "e  + V'  ~e 
i it  t in  n 


(A-13) 


For  the  smooth  sphere  model, 

\ 2 


it 


-el 


in 


cos  x = 


(A-14) 
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In  terms  of  the  collision  parameter  (b)  and  collision  cross-section  (a) 
(Fig.  A-2), 


cos  x 


2 2 
sin  a - e cos  6 

77TT7  2 2nv  1/2  ’ 

(sin  0 + e cos  0) 


where 


Hence , 


and 


sin  0 = b , 
b = b/0 


cos  x = 


v2 


(l+e)b2  -e 


r/1  2,-2  , 2,1/2  * 

[ (1-e  )b  + e ] 


= (l-e2)b2  + e2 


For  e=l  (adiabatic) 


-2 


cos  X = 2b  - 1,  , b £ 1, 
■ 1 , b > 1, 


V ■ V 
i i 


For  e=0  (inelastic) 


cos  x = b , b <_  1 , 

= 1 , b > 1 


vi  - bvi  • 


(A-15) 


(A-16) 


(A-17) 


(A-18) 


(A-19) 


Generalizing  the  smooth  sphere  model,  for  adiabatic  collisions  with 

11 

spherically  symmetrical  interaction  potentials,  Hirschf elder , et_  al^.  gives 
for  the  CM  scattering  angle, 
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00 


where 


X(b,g) 


TT-2b 
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m 


dr 
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(1 


± 


1 2 
2^g 


Tfi 


(A-20) 


r = minimum  distance  between  the  molecules  during  the  collision 
m 

4>(r)  = interaction  potential 

g = initial  relative  speed. 

The  elastic  sphere  (e=l)  is  then  given  by, 

<t>(r)  = 00  , r < a , 

- 0 , r > 0 , 


r = O , b < O , 
m — 

= b , b > a . (A-21) 

and  yields  Eqs.  (A-18)  upon  integration. 

In  summary,  for  adiabatic  collisions  there  is  no  energy  loss  (by 
definition),  the  speeds  in  the  CM  system  are  unchanged,  the  collision  is 
completely  described  by  Eq.  (A-20)  for  spherical  potentials,  X is  the  same  for 
both  particles  and  the  impact  parameters  are  unchanged.  Non-adiabatic 
collisions  have  been  modeled  by  smooth  spheres  with  the  coefficient  of 
restitution. 

LES-8/9  Scattering  in  the  Laboratory  System:  The  speed  of  the  ambient 

molecules  was  taken  to  be  equal  to  the  speed  of  the  launch  vehicle.  The  speed 
of  the  desorbed  molecules  with  respect  to  the  emitting  surface  was  taken  as 
the  mean  thermal  speed  evaluated  at  the  surface  temperature  (100°C) . Table  2 
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shows  the  desorbed/ambient  momentum  ratios  for  the  LES-8/9  launch  at  an  altitude 
of  132.55  km  (305  seconds  after  ignition)  for  a range  of  desorbed  molecular 
weights.  While  the  desorbed  molecules  are  more  massive  they  travel  sufficiently 
slow  so  that  the  ambient  molecules  have  more  momentum. 


TABLE  2 

MOMENTUM  OF  DESORBED  AND  AMBIENT  MOLECULES 
(LES-8/9  launch  - 132.55  km) 


Desorbed 
Molecular  Weight 

1000 

500 

100 

50 


Desorbed 

Thermal  Velocity* 

_2 

8.8866  x 10  km/sec 
1.255  jc  10_1 
2.8101  x 10-1 
3.9742  x 10_1 


Desorbed /Ambient 
Momentum  Ratio** 

0.6149 

0.4348 

0.1944 

0.1375 


T = 100°C 

W = 27.5  , V = 5.255  km/sec 
amb  amb 


Figure  11-3  schematically  shows  the  scattering  of  a desorbed  molecule 
emitted  at  45°  to  the  surface  normal  by  an  ambient  molecule  traveling  parallel 
to  the  surface.  The  figure  corresponds  to  a typical  collision  during  the 
LES-8/9  launch  when  both  molecules  are  modeled  by  perfectly  elastic,  smooth 
spheres.  As  shown  in  the  preceding  subsections  for  this  model,  the  desorbed 
molecules  are  scattered  uniformly  over  a sphere  centered  at  v^  + vCM,  i.e.. 
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isotropically  with  respect  to  the  center  of  mass  system.  Figure  II-3  shows 
that  the  assumption  of  isotropic  scattering  in  the  laboratory  system  will 
underestimate  the  number  of  molecules  scattered  in  the  forward  and  downstream 
directions  while  over-estimating  the  number  scattered  in  the  backward  and  up- 
stream directions.  The  equations  developed  in  Sections  II  through  IV  therefore 
overbound  the  molecular  irradiance  on  upstream  surfaces  of  the  spacecraft. 
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APPENDIX  B 


EVALUATION  OF  SOME  INTEGRALS 

In  Section  III-A  (p.  13)  it  was  shown  that  the  molecular  irradiance  due  to 
backscattering  from  a cylindrical  beam  was  a function  of  the  following  integral, 


I^e)  = / dx  exp(-x)  |l  — 
n L 


, 2 2.1/2 

(x  +e  ) 


It  is  desired  to  evaluate  Eq.  (B-l)  for  e <<  1.  From  Gradshteyn  and  Ryshik 
Eq.  (B-l)  is 


(B-l) 
5 


\(e) 


1 - 


j e^H-^e)  - N1(e)  - ej  | 


(B-2) 


where  = Struve  function  first  order 

N^  = Neumann  function  first  order 


By  definition 


Hx(e) 


00  (-)  Or) 

E 


2+2m 


m=o 


r (m+3/2)  r (mf5/2) 


ttN1(e) 


2J1(e)  (£n  | + C)  - | - | 


, \k/£\ 
00  (-)  (-?) 

E 


l+2k 


k=l  k!  (k+1) 


1+k  k .. 

E — + F — 
m=l  m m=1  m 


(B-3) 


where  C = Euler's  constant  and  is  the  Bessel  function  of  the  first  kind. 


Writing  out  the  first  few  terms  of  (B-3)  and  substituting  in  (B-2)  gives 
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!]_(£) 


= 1 + e — 


IT 


(~)2 
_ T_ 

r(|)  r(|) 


p 4 

(§r 

r(|)  r<|) 


+ . 


-f  Jl(E)  Mnf  +C,  +i  (f)  +i  (f) 


h <!>3  <f>  +' 


For  e « 1 


I^e)  - 1 + e -j  e (±  |)  +0(e2) 

or 

I1(e)  « e 


(B-4) 


(B-5) 


(B-6) 


For  £ small  I^(e)  is  determined  mainly  by  values  of  the  integrand  near  x = 0. 
Hence,  bounding  exp  (-x)  by  unity  and  integrating  gives  the  first  term  in  the 
expansion  of  the  exact  solution. 

In  Section  IV-A  (p.  20)  it  was  shown  that  the  far-field  solution  for  the 
molecular  irradiance  was  a function  of  the  following  integral, 


I 

03 


(e  • e + e * e ) 
r r j « s r 03  r 

II  do)  cos0  

emitting  6 (1  + eU)  eg) 

cone 


(B-7) 


Using  the  coordinate  system  shown  in  Fig.  IV-1  it  can  be  shown  that 


e • e = cosy 
s r s 


1 + e 

03 


e 

s 


a + b cos^ 


e “ e = c + d cos d>  -I-  e sind 
03  r e e 


(B-8) 
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where 


a = 1 4*  cos0  cos0 
s e 


b =*  -sin0  sin0 
s e 


c = COS0  COS0 

e 


d = sin0  cos <b  sin0 

e 


e =*  sin©  sin (b  sin0 

e 


Eq.  (B-7)  becomes 

0 

I * / d0  sin0  cos0  1.(0  ) 
0)  J e e e <p  e 

o 

where 


I A6J  * I-  + I9  + I 

<p  e 12  3 

2tt  d<j> 

T1  ■ <cos*s  + c>  / a + /cos* 
o e 

2tt  d<£  coS(j) 

I9  = d / | , — 

2 JQ  a + b cos^ 

2tt  d(J)  sin<j) 

I-  = e / — f— T 

3 J a + b cos(J) 
o 

Evaluating  1^,  and  gives 


1 ee> 


2tt(cosy  + c)  , 2tt  — 

S _l  2 TT  — - 

9 9 1/2  + z b „ , 1/2 

(aZ  - bZ)  (a^  - b^) 


Substituting  (B-10)  in  (B-9)  gives 


(B-9) 


(B-10) 
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I = 2tt 
0) 


sinQ  cos(ft  sin  6 

sin0  2 

s 


( . sin0  cosd) 

+ cosy  + — — x 

\ s sin0 


)[- 


/ 1 + COS0  V. 

( —) 

' COS0  + COS0  / 
s J 


cos0  — cos0  £n 
s 


+ ( cos0  + 


sin0  coscf)  COS0^ 
sin0 


sin20  . 2, 

— — — — cos0s  (1  — COS0)  + cos 


/ 1 + cos0  \1  ) 

9 £n  ( k— V-J  (-(E-ll) 

s \cos0  + cos0  / J 


Letting 


I = 7 f 
W I 9 


(B-12) 


where  for  a Lambertian  surface 


N = tt!  sin  0. 


f = - 


sin0  cosft 
sin0 


0/  , sin0  cosd>\ 

+ 2 (co s1rs  + sine  ' V 

S 

sinS-  COS(J)  COS0 
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1 — COS0  COS  0 g ( s 
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1 ~ cos0 

sin20 
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2 cos  0 

s 

sin^0 


/I  + COS0s  V | 

\ COS0  + COS0  / j 
S J ' 


£n 


(B-13) 


Multiplying  out  and  collecting  terms  simplifies  (B-13)  to 


f = 


COS0  + cosy 

' c 

1 + COS0 


(B-14) 
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NOMENCLATURE 


A 

b 

C 

C 

-> 

e 

H 

I 

K(a,g) 

K 

n 

lyl 

o 

M,m 

N 

N 

o 

r 

s 

S 

t 

v 

V 

V 
W 

V 


amb 


Area 

{ Surface  dimension  (Fig.  IV-4) 
\ Collision  parameter 


cm 


cm 

cm 


Surface  contamination 
Time  rate  of  change  of  C 
Area  normal  unit  vectors 
Molecular  irradiance 
Molecular  specific  intensity 
Geometric  factor  for  scattering  from  a plane 
Knud sen  number 
Lengths  (Fig.  IV-4) 

Molecular  mass 
Emitted  molecular  flux 
Avogadro T s number 
Radial  coordinate 


-2 

gm  cm 

-2  -1 
gm  cm  sec 


molecules  cm  ^ sec  ^ 
molecular  cm  ^ sec  1 str  1 


Distance  between  emitting  and  receiving 
areas 

Nondimensional  s 
Time 

Velocity 

Volume 

Velocity  of  satellite 
Molecular  weight 

Time  rate  of  changes  of  paint  coating 
mass  fraction 


cm 


gm 


molecules  cm  ^ sec  ^ 


molecules  moles 


-1 


cm 

cm 

cm 

sec 


cm  sec 
3 

cm 
km  sec 


-1 


-1 


gm  mole 
-1 


-1 


sec 


x 


Nondimensional  r 
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O?  "CD  P 


y Nondimensional  r 

r 


Nondimensional  lengths  (Eq.  49) 

Paint  thickness  mils 

G Eq.  (17) 

Y Eq.  (24) 

X Mean  free  path  cm 

! Paint  density  gm  cm 

Eq.  (16)  cm 

0 Polar  angle  deg. 

4>  Azimuth  angle  deg. 

X,Xf  Center  of  mass  and  laboratory  scattering  deg. 

angles 

u)  Emitting  solid  angle  str 

u)f  Receiving  area  solid  angle  str 


Subscripts 

e Emitting  area 

r Receiving  area 

rel  Relative 

amb  Ambient 
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